Stable "antiferromagnetic" vortex lattice imprinted into a type-II superconductor 
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In type-II superconductors, where vortices and antivortices tend to annihilate, only a "ferromag- 
netic" vortex lattice, with the same orientation of vortex magnetic moments, is usually formed in 
a homogeneous external magnetic field. Using the time-dependent Ginzburg-Landau formalism, 
we demonstrate that a checkerboard vortex-antivortex lattice ("antiferromagnetic vortex lattice"), 
imprinted onto a superconducting film by a periodic array of underlying clockwise and counterclock- 
wise microcoils generating spatially periodic positive and negative magnetic field pulses and then 
trapped by an array of artificial pinning centers, remains stable even after the imprinting magnetic 
field pulse is switched off. 



Magnetic moments in solids can order ferro- or antifcr- 
romagnetically In a certain analogy with that, in 

type-II superconductors the ordered vortex state can be 
built up from entities of the same or the opposite polarity 
- vortices and antivortices. Such vortex lattices can be 
referred to as "ferromagnetic" and "antiferromagnetic" 
vortex states. In this terminology, the Abrikosov vor- 
tex lattice [i| is an example of a "ferromagnetic" state. 
During the last decade, vortex patterns consisting of su- 
perconducting vortices and antivortices have attracted 
significant theoretical and experimental efforts. A num- 
ber of stable vortex-antivortex (VAV) configurations have 
been predicted and experimentally detected in symmet- 
ric mesoscopic superconducting samples subjected to a 
homogeneous magnetic field [5|-|7| . VAV pairs can be nat- 
urally created and stabilized by spatially inhomogeneous 
fields of micromagnets jl - 11 1 , leading to various commen- 
surability effects and novel stable VAV configurations in 
superconductor- ferromagnet hybrids with regular arrays 
of magnetic dipoles (il. [l2 . Il3 1 . 

Whereas VAV lattices were induced using an array of 
micromagnets, the creation of a VAV state that remains 
stable without permanent presence of external spatially- 
alternating magnetic field is still an open problem. An 
appealing idea relying on superconducting relaxation dy- 



namics was formulated in Rcfs. [1J, |15|. This pro- 
posal considers vortices and antivortices that are sponta- 
neously generated during the recovery of superconductiv- 
ity after a sample was heated locally with a laser pulse. 
In the presence of strong pinning centers some of those 
vortices and antivortices can be frozen and persist prac- 
tically indefinitely. In this context it was suggested [16 1 



that an analysis of the superconducting phase recovery 
accompanied by the creation of vortices and antivortices, 
provides a convenient tool to test the Kibble-Zurek cos- 
mological scenario 17, [Hf for nucleation of topological 



defects during a symmetry-breaking phase transition. 

In this Letter, we consider the use of a specially de- 
signed pulse of spatially periodic magnetic field for im- 
printing a regular VAV lattice, which, as we predict, re- 



mains stable after the removal of the external magnetic 
field used for VAV imprinting. In particular, we will 
demonstrate that due to mutual cancelation of forces, 
which act on a vortex (antivortex) in the equilibrium 
state of such a lattice, the antiferromagnetic vortex state 
can be stabilized even in the case when pinning centers 
in the superconductor are relatively weak. 

To describe the vortex dynamics, we apply the time- 
dependent Ginz burg -Landau (TDGL) formalism as pre- 
sented in Rcfs. |l9l lioj. For the sake of convenience, 
we use dimensionless variables by expressing lengths in 
units of \/2£ and currents in units of $o£/ (v^t^oA 2 ), 
where $o — wh/e is the magnetic flux quantum, fio is 
the vacuum permeability, £ is the coherence length and 
A is the penetration depth. The unit of magnetic field, 
<I?o/(47r£ 2 ), is half the second critical field. Our unit of 
time is 2wtgl, where tql is the Ginzburg-Landau relax- 
ation time and u — 7r 4 /[14£(3)] with the Riemann 
zeta function. 

We consider a thin infinite superconductor film sub- 
jected to a pulse of spatially periodic negative and posi- 
tive magnetic fields [Fig. 1(a)]. In the model, this field is 
generated by an underlying square array of circular quasi- 
one-dimensional current loops with diameter D, located 
at a distance d\ below the middle of the superconductor 
layer. The current direction changes from one loop to 
another in the checkerboard order as shown in Fig. 1(a), 
thus generating a periodic array of local magnetic fields of 
the opposite polarities. It is assumed that due to an ap- 
propriate arrangement of the leads connecting the loops 
to each other, the magnetic field induced by these leads 
is negligible. In the present calculations, the spatial pe- 
riod of the structure is L — 15, the superconductor film 
thickness is d = 0.3 and the Ginzburg-Landau parameter 
of the superconductor is k = A/£ = 1. The current pulse 
in the loops has trapezoidal shape: the current I(t) lin- 
early increases from to / max in the time interval from 
t = to t = ti, remains constant I(t) = I ma x in the 
time interval fa, £2) and then linearly decreases down 
to within the interval fa, £3). A periodic lattice of 
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up and down out-of-planc magnetic fields, induced by 
such a current pulse, creates in the superconducting layer 
a square lattice of vortices and antivortices located just 
above the loops with counterclockwise and clockwise cur- 
rent directions, respectively. Vortices and antivortices, 
arranged to form a periodic checkerboard-like square lat- 
tice, remain in equilibrium even after completely switch- 
ing off the applied magnetic field. This equilibrium is 
obviously unstable in an ideally uniform superconduct- 
ing film: an arbitrarily weak perturbation would destroy 
the lattice due to the vortex-antivortex attraction and 
subsequent recombination of the VAV pairs. However, 
in reality superconducting films often contain numerous 
naturally occurring pinning centers, which may prevent 
such a collapse of the VAV lattice. Alternatively, in a 
more controlled way, a periodic array of artificial pinning 
centers can be intentionally induced, e.g., in the form of 
a periodic array of "blind" or complete antidots [U [22| . 

Here we consider square array of artificial pinning cen- 
ters formed by indentations of depth 6 -C d made in the 
superconducting film of thickness d (blind antidots), as 
shown in Fig. 1(a). The number and periodicity of the 
pinning sites corresponds to the number and periodicity 
of the current loops, but the arrays of loops and blind 
antidots are shifted with respect to each other. We have 
analyzed shifts along the side of a unit cell, over a dis- 
tance denoted by s x , and shifts along the diagonal of the 
unit cell, over a distance Sd- 

Figures 1(b) and 1(c) illustrate the formation of a VAV 
lattice in the case of 6 = 0.06, di = 0.5, D = 1, s x = 2, 
-fmax = 4, ti = 10, t 2 = 210, and t 3 = 250. Ramping 
up the periodical array of up and down magnetic fields 
during the time interval (0, t±) strongly suppresses the 
order parameter \ip\ just above the centers of the current 
loops. Due to the radial gradient of the magnetic field, 
induced by a current loop, a VAV pair is generated at 
each loop. This process is quite similar to the VAV pai r 
generation by a magnetic dipole considered in Ref. [231 ] . 
As the current pulse reaches its plateau, some recom- 
binations take place: the vortex induced by the field 
of a current loop with counterclockwise current remains 
pinned to this loop, while the corresponding antivortex 
is driven away from the loop and recombines with a vor- 
tex pushed away from one of the neighboring loops with 
clockwise current. After that, single vortices and antivor- 
tices, which correspond to one flux quantum <£>o each, 
are left at the sites of the counterclockwise and clockwise 
current loops, respectively. Due to the aforedescribed 
processes, a VAV lattice, pinned to the current-loop ar- 
ray, is imprinted into the superconductor [see Fig. 1(b)]. 
This lattice pertains till the end (t = t^) of the plateau in 
the applied field pulse. During the time interval fa, £3), 
when the applied magnetic field decreases and eventually 
vanishes, the vortices and antivortices gradually drift to- 
wards the nearest pinning sites, resulting in the formation 




FIG. 1. (a) Unit cell of a periodic structure, which contains a 
superconducting film and a square array of f2-shaped current 
loops of diameter D at a distance d\ from the film. When 
applying a current I(t) as shown with green arrows, the mag- 
netic field induced by the loops alternates in checkerboard 
order. Spatial modulations of the superconductor thickness 
provide an array of vortex-pinning centers, which are shifted 
along the :r-axis by a distance s x with respect to the cen- 
ters of the nearest current loops, (b) At the end of the pulse 
plateau a lattice of vortices and antivortices, seen as local 
suppressions in the square modulus of the order parameter 
\ip\ 2 , is created by the magnetic field B z of the loops. The red 
(blue) arrows indicate the direction of the magnetic moment 
of vortices (antivortices). (c) When the applied current de- 
creases and eventually vanishes, the vortices and antivortices 
drift towards the nearest pinning centers. This results in a 
stable VAV lattice characterized by an "antiferromagnetic" 
distribution of the induced magnetic field B z . 



of an "antifcrromagnctic" vortex state, which is shown in 
Fig. 1(c) and which remains stable for an arbitrarily long 
time after the applied field has been switched off. 

Figure 2 summarizes the results of the TDGL simula- 
tions for the structure shown in Fig. 1(a) at S = 0.06, 
di = 0.5, D = 1, , ii = 10, t 2 = 210, and i 3 = 250. 
Depending on the magnitude of the current in the loops, 
/max, and the horizontal (s x ) or diagonal (sd) shift be- 
tween the centers of the current loops and the pinning 
centers, different vortex configurations, shown in Fig. 2 
with different symbols, appear to be stable after switch- 
ing off the magnetic field pulse. In the case of zero shifts 
s x and Sd, no VAV lattice is formed in the film for cur- 
rent magnitudes lower than 7 max ~ 3.6, even though at 
J mas < 3.6 the applied magnetic field pulse can lead to 
strong local suppressions of the order parameter. Inter- 
estingly, a moderate nonzero shift s x facilitates nucle- 
ation and spatial separation of vortices and antivortices 
from these suppressions. As follows from Fig. 2, this 
can result in the formation of a stable VAV lattice even 
for current magnitudes I max slightly below 3.6. A sim- 
ilar, but somewhat less pronounced effect, is seen to be 
caused by a diagonal shift Sd- 

At 3.6 < / max < 4.6, an antiferromagnetic lattice is 
stabilized both at zero shifts s x or Sd and also at mod- 
erate nonzero values of these shifts. Our calculations 
show that for this range of / max the vortex pattern cre- 
ated at the end of the pulse plateau (t = t 2 ) has exactly 
one single flux quantum vortex or antivortex near the 
center of each current loop, thus forming a VAV checker- 
board lattice. Taking as a unit cell the square with side 
L [see Fig. 1(a)], this leads to two VAV pairs per unit 
cell. At moderate shifts s x or Sd, when switching off the 
applied magnetic field, this lattice simply shifts to the 
pinning-site array. However, when the distance between 
the centers of the loops and the pinning sites exceeds a 
critical value, complete VAV annihilation, rather than 
the antiferromagnetic-lattice stabilization, occurs after 
switching off the applied magnetic field. 

At larger I max , i.e. for higher magnetic fields, each 
current loop generates more than one VAV pairs due to 
the application of the magnetic field pulse. As a re- 
sult, a giant vortex (antivortex), which corresponds to 
two or more flux quanta 0, HH, or a dense cluster 
of single flux quantum vortices (antivortices) is accom- 
modated by most current loops at the end of the pulse 
plateau. Ramping down the applied field during the time 
interval (t 2 , £3) leads to dissociation of these giant vor- 
tices (antivortices) or vortex (antivortex) clusters into 
smglc-<I>o vortices (antivortices). The further evolution 
of the system is determined by an intricate interplay of 
the vortex- (anti) vortex interactions as well as the inter- 
action of vortices and antivortices with the pinning po- 
tentials and the gradually decreasing applied inhomoge- 
neous magnetic field. This leads to a rather complicated 
recombination dynamics of the vortices and antivortices. 
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FIG. 2. Results for stabilization of a VAV lattice in the struc- 
ture shown in Fig. 1(a) for different current magnitudes / ma x 
and shifts s x (a) or Sd (b) between the centers of the cur- 
rent loops and the pinning centers. The dotted lines approxi- 
mately indicate the boundaries between the ranges of param- 
eters which correspond to stabilization of a VAV lattice with 
two (V = 2, light blue filled circles) or one (V = 1, dark blue 
semi-filled circles) VAV pairs per unit cell, and the ranges of 
parameters which correspond to eventual annihilation of pairs 
(V = 0, empty circles). 



After that, generally one of two outcomes is realized: ei- 
ther we again obtain a stable VAV checkerboard pattern 
with two VAV pairs per unit cell as before, or full re- 
combination has taken place. Remarkably, a moderate 
diagonal shift Sd of the current loops with respect to the 
pinning sites can promote the antiferromagnetic- vortex- 
lattice stabilization for this range of I max . No such effect 
is found for a shift s x along the side of the unit cell. The 
following qualitative difference between diagonal and hor- 
izontal shifts seems to be important: in the former case, 
the pinning sites lie on the lines that connect loops with 
the same direction of current, while in the latter case 
the pinning cites arc shifted from a current loop towards 
a loop with the opposite current direction. As further 
seen from Fig. 2, besides the two "limiting cases" (sta- 
bilization of a lattice with two VAV pairs per unit cell 
or full annihilation of pairs), sometimes an "intermedi- 
ate" configuration is observed for parameter values near 
the boundary between the regions corresponding to the 
aforementioned two "limits" . This configuration has one 
VAV pair per unit cell, and an example of this pattern 
will be given below. 

Large shifts between the pinning lattice and the cur- 
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rent loop lattice impede the stabilization of the VAV 
lattice. In order to avoid this, the density of artificial 
pinning sites can be increased so that a pinning site is al- 
ways nearby for the vortices generated by a current loop. 
An example of such a structure is shown in Fig. 3(a). 
As in the cases analyzed above, stable antiferromagnetic 
vortex lattices in this structure may contain either two 
VAV pairs per unit cell or one VAV pair per unit cell. 
The latter is illustrated by Fig. 3b, showing the evo- 
lution of the magnetic field patterns, corresponding to 
vortices and antivortices, for the particular set of param- 
eters, S = 0.09, d x = 1, D = 2, I max = 5.5, h = 10, 
t-2 = 160, and ta = 260. By the end of the magnetic-pulse 
plateau (t = 160), giant vortices (antivortices), carrying 
two flux quanta, are formed at each current loop. When 
decreasing the applied magnetic field, these giant vortices 
and antivortices still persist up to t rj 244 [see Fig. 3(b)]. 
However, with further reduction of the applied field (at 
t rj 252), they dissociate into single-$o vortices and an- 
tivortices, which then partially recombine [see the snap- 
shots for t = 254 and 255 in Fig. 3(b)]. Since the pinning 
potentials in the structures under consideration are rela- 
tively weak, they do not prevent the corresponding mo- 
tion of vortices and antivortices. This means that these 
pinning centers can provide fine tuning in stabilizing a 
state with coexisting vortices and antivortices. The cru- 
cial prerequisite for such a stabilization is mutual cance- 
lation of vortex- vortex and vortex-antivortex interactions 
in a regular antiferromagnetic VAV lattice. In the case 
under consideration, the order parameter finally settles 
into a stable state with one VAV pair per unit cell, cor- 
responding to periodic VAV chains, separated from each 
other by a distance L [see the corresponding magnetic 
field distribution at t = 350 in Fig. 3(b)]). 

Each current loop, independently of the sense of the 
current in it, can generate both vortices and antivortices 
in the superconductor film. This suggests that a stable 
VAV lattice can be formed also in a periodic structure 
where the current direction is the same for all the current 
loops. The results of our calculations for those structures 
(see supporting information) confirm such a possibility. 
Moreover, these results demonstrate that (at least for 
certain sets of the relevant parameters) stabilization of 
a VAV lattice can be possible even when the number of 
vortices and antivortices, induced by each current loop in 
the course of the magnetic field pulse, is much larger than 
one. At the same time, we find that the most promising 
and reliable regime for the experimental realization of 
the imprinted antiferromagnetic vortexantivortex lattices 
that remain stable after the external magnetic field is 
switched off, corresponds to the generation of just one 
VAV pair by each current loop. 

Since the unit of length, used in our calculations, is pro- 
portional to the coherence length £, the examples given 
above can correspond to different geometric parameters, 
depending on a particular superconducting material and 




FIG. 3. (a) Unit cell with four current loops and a relatively 
high density of pinning centers in the superconductor film, 
(b) Snapshots of the magnetic field B\ z induced by the su- 
perconducting film at different time moments t. 



temperature. Thus, for the coherence length £ rj 250 nm 
found in Ref. [H for an Al film at T = 0.63T C , the 
structure shown in Fig. [3] corresponds to the unit-cell 
size about 5.3 /zm, the superconducting-layer thickness 
100 nm, the diameter of current loops 700 nm and the 
spacing between the loops and the superconducting layer 
about 350 nm. For the depth and lateral size of the inden- 
tations, which serve as pinning sites, we then obtain the 
values about 30 nm and 600 nm, respectively. Of course, 
all the indicated values would appear larger (smaller) for 
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higher (lower) temperatures. 

Experimentally, stable imprinted VAV lattices, the- 
oretically predicted in this paper, can be directly vi- 
sualized in nanopatterned systems (similar to the one 
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or 



shown in Fig. [TJ by using scanning SQUID 
Hall probe [30( microscopies. 

To conclude, we have proposed the application of a 
pulse of spatially periodic up and down magnetic fields 
to imprint a vortex-antivortex lattice on a superconduct- 
ing film. The stabilization of the imprinted VAV pattern 
can be provided by an artificial array of pinning centers, 
so that the antiferromagnctic vortex pattern remains sta- 
ble in zero magnetic field, after the magnetic field pulse 
is turned off. The presence of pinning centers (especially, 
of their periodic arrays) in the superconducting layer is 
a necessary but - in general - not sufficient condition 
for the antiferromagnetic-vortex-lattice stabilization: it 
is sensitive also to the magnitude of the applied mag- 
netic field pulse as well as to its spatial and temporal 
shape. Nevertheless, our results clearly demonstrate that 
the relevant parameters, which assure creation of stable 
antiferromagnctic vortex lattices with a relatively high 
density of vortices and antivortices, can vary in a rather 
broad range. Moreover, the stabilization of the antifer- 
romagnctic vortex lattices is even easier in presence of a 
denser periodic array of artificial pinning centers. 
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